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Abstract 

The flow of granular material through an orifice is studied experimentally as a function of force 
F pushing the flow. It is found that the flow rate increases linearly with F — a new, unexpected 
result that is in contrast to the usual view that F, completely screened by an arch formed around 
the orifice, has no way of altering the rate. Employing energy balance, we show that this behavior 
results mainly from dissipation in the granular material. 
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Although granular materials such as sand, rice etc. are familiar in daily life and industrial 

handling, their static and dynamic behavior is only beginning to be understood l| . In 

these systems, arch formation is an important effect, which may remarkbly influence the 

mechanical properties of these materials^]. A classic example of a dynamic arch can be 

found in the simple system of a silo (fig. la), in which mass is discharged from an orifice. 

Since ancient times, it has been supposed (and used for making "hour glasses") that the 

discharge rate Q of the granular orifice flow (GOF) is independent of the height H of the 

silo column — a phenomenon commonly believed to be due to the so called "free fall arch" 

fl 

(FFA) formed around the orifice[3(]. Moreover, if the orifice diameter D decreases, the FFA 
may jam into an static arch, at which point the flow stops. The jamming transition was 



studied in J4] for the 2D case. Specifically the FFA is described as a semispherical surface 
spanning the orifice, the total energy of which is taken as minimal in the "minimum energy 
theory" jj, and the normal stress taken as vanishing in the "hour glass theory" jf]]. 

These theories represent the most advanced continuous modeling of GOF, from which 
the empirical Beverloo formula for Q can be derived^. The FFA assumed here is rather 
unusual. It can resist all force F given by the weight of the column H, and any load L on 
its top, such that Q is fully protected from being influenced. Clearly, this basic assumption 
can be directly checked, by measuring the rate Q as a function of F. We performed this 
measurements and found a linear dependence (if F is not too small), 

Q = A{l + aF) (1) 

instead dQ/dF = 0. This demonstrates that the properties claimed for FFA in 

BB 

are, 

at the least, not rigorously valid. The difference between (CQ) and the Torricelli's law for 
ideal liquids: Q ~ \fF is clearly a result of granular arch and dissipation. Note that to 
measure the variation of Q with F, we need to remove the Janssen effect 
L and the weight of sand are typically redirected onto the side wall of silo, steering F to a 
saturation value F sat determined by the frictional property between sand and wall, including 
the history of the preparation. The screening of F can be removed, for instance, by allowing 
the side wall to move freely and vertically (see fig. lb), then F is a sum of the load L, the 
weight of sand M, and the wall My/ - a quantity that can be experimentally controlled and 
varied. 



3, |8j , that the load 



The variation of Q with the orifice diameter D has been studied for a long time 
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The Beverloo formula, 



Q = Cp^g(D-D a 



,5/2 



(2) 



though well known, needs as a topic of current interests, validation and a clarification of 
its limitation, see e.g. [10;]. Here, g is the gravitational acceleration, C a parameter, and 
D a a correction to the orifice diameter of an empty annulus, necessary because the flow has 
a tendency to avoid the periphery of the orifice S] ■ More recently, Sheldon and Durian 
measured the variation of Q with the orifice tilt angle 9 11] . Neither dependence (not to 
mention their combination) is well understood, posing a serious challenge to any granular 
theory. We shall consider only horizontal orifice planes (6 = 0), focusing our attention on 
the variation of Q with the pushing force F. To our best knowledge, this property has not 
been addressed before. There is a simulation which shows an increasing of the GOF velocity 
with column height 12], which seems similar to (HI), but its authors attributed this effect 
to the choice made in the simulation of the damping parameter, and did not take it as a 
general behavior of GOF. Moreover, strong vibration may break up the arches, dramatically 
altering the flow behavior, even rendering it to become the Torricelli one [3, [l^.Also, for 
dilute granular flows studied in a deviation from (CD) is expected, because no arches are 
formed. 

Experimental setup is shown in fig. lb, where the side wall is a PVC cylinder of inner 





FIG. 1: Granular flow from a bottom orifice with (a): fixed wall, and (b): vertically movable wall. 
(See text). 
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FIG. 2: g 2 / 5 versus D. Full straight is the best linear fit. Inset shows a measured discharging 



M(t). 



diameter <\> = 10.4 cm, weight Myy = 1.16 kg. It is horizontally restricted by pulley wheels 
such that it can be moved freely in the vertical direction. The granular material consists of 
nearly mono-disperse glass beads of the diameter d ~ 3 mm. Bottom orifice is of D = 14, 
16, 18, 20, 22 mm. Load L ranges from to about 40 kg. For each experiment, the flow 
starts with a granular column of the height Hi = 51 cm, prepared by direct pouring, and 
ends with Hf = 37 cm, the total flowed mass is about Mo = 1.8 kg. The bulk density p of 
the granular sample is about 1.5 g/cm 3 . GOF rate Q is measured by continuing recording 
(15 records/s) the flowed mass M = f* Qdt' with an electronic balance of the precision ±0.5 
g. It is observed that M (t) is a fairly good straight line, indicating a constant rate Q given 
by its slope (see e.g. fig.2 inset). For every measurement the flow rate is found reproducible, 
with a fluctuation less than 1%. First, we performed measurements with a fixed cylinder 
without load (L = 0), and obtained a rate dependence W au(-^) that f°ll° ws the Beverloo 
law (J2J) excellently, as shown in fig.2. Best fitting gives D a ~ 3.7 mm and C ~ 0.51, close 
to the values 4.5 and 0.58 given in the book j^]. 

Next, we release the cylinder, allowing it to move vertically, such that the silo screening 
effect is eliminated. In this situation, the bottom force F pushing the flow is given by the 
total mass of the setup, which varies slightly with the time, M = M (t). During the flow, 
which starts at t — and terminates at tf, the temporally averaged pushing force is 




Fdt = L + M w + 



(Hi + Hf) 
2 (Hi - H f ) 



(3) 
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FIG. 3: Measured GOF rate Q versus average pushing force (F) (symbols), and their best linear 
fit (straights), for various orifice diameters D. 

Because no noticeable influences of this slight temporal dependence of F on the rate Q has 
been observed, (the measured J Qdt is always a good straight line as given in fig. 2 inset), we 
shall assume F ~ (F) in what follows. Fig.3 shows the measured variation of Q with (F), 
all straight lines, confirming Eq.([T]) with a > 0. Note that although the load, the cylinder 
and the upper part of the granular column go down slowly during the flow as a whole, no 
relative motion among them is observed. 

To check the validity of the Beverloo law (j2J) for pushed GOF, we plot Q 2 ^ versus D for 
various forces (F) in fig. 4a. Again linearity was obtained, implying a good agreement. By 
best linear fit, we found that both Beverloo parameters C, D a increase with pushing force 
(fig.4b and c), with the former showing a linear behavior, 

C = Co (1 + (F) /F ) (4) 

(where Cq = 0.44 and F = 35 kg), while the latter is nonlinear, displaying a saturation 
behavior as (F) increases, approximately as 

n - n - D °° ~ — 
a 00 1 + (F) /F* [ ' 

(where is its biggest value at the large force limit (F) — > oo, and Dq the smallest one 
at the zero force limit). Best fit yields F * ~ Fq, ~ 8.3 mm, Dq ~ 3.2 mm, see fig. 4c. 

In the Beverloo law, D a accounts for clogging, which happens when the orifice is meso- 
scopic, of a size no more than a few grains. Our measurements show that clogging is slightly 
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FIG. 4: (a): Q 2 / 5 versus D for various pushing force (F). Symbols are the data in fig. 3, and full 
straights their best linear fits. Their Beverloo parameters are plotted in (b) and (c) (symbols) as 
functions of (F), where full curves are the model (|4l5p . 



more likely with increased force. Note that as a finite size effect clogging cannot exist in 
the macroscopic limit, D » d, where the Coulomb yield law holds. This means D a must 
saturate at a mesoscopic size. 

In the above generalization of the Beverloo law which includes the pushing force, varia- 
tions of C, D a as modeled by (Hll5|) is considered. Unfortunately, this way of generalization 
obscures the linear relation between Q and F. Alternatively, we may consider a different 
generalization that keeps the linearity: Letting A and a in (pQ) vary with D as 

A = p^C (D-D f\ (6) 



When the model constants C , D , F take the same values as in (j4f5]) . the two models 
become identical in the limit (F) — >■ or D » D oo l . As the difference between them 
is very small (see fig. 5a), it can not be decided with the present data which version of 
generalization is more appropriate. Using the values for A, a given by the straight lines of 
fig. 3, we obtain a good agreement with f[rJfT|) . if D\ ~ 13 mm (fig. 5b, c). 

It is worth noting that the above empirical formulas for Q (F, D) are invalid in the limits 
of small pushing forces F — > and clogging orifice diameter D — > D c \ og > D^, where Q is 
vanishing but the formulas do not show it. The diameter D c i og , at which GOF is clogged, 
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FIG. 5: (a): Comparison between the model (|4|5jl (circls) and the model (|6|7p (curve). Inset is an 
ampified figure for small D. (b): variation of a , (c): of A, with L>. Symbols are measurements, 
and curves are the model (|6|7|) . 



is F dependent, and not a parameter of these models. For the experiments of this work, we 
have D\ > D c \ og > Doo > D . That -D c i og could not be obtained with the models is probably 
due to the fact that the transition from flow to clogging is not a sharp one, but occurs via 
an intermittent flow. The above values for -Do,oo is well below the transition at which no 
flow is possible, while Di belongs to the intermittent regime, which starts from D ~ 12 mm 
for L = 0, and ~ 14 mm for L = 40 kg. Moreover, in addition to causing the flow, F also 
increases clogging, as it jams the grains around the orifice and FFA. At D ~ Di, these two 
effects compensate each other, and we have a ~ 0, implying Q's independence of F. 

The dependence of the GOF rate on the pushing force F can be qualitatively understood 
employing the energy balance: The power Wi n injected into the granular system is equal to 
the sum of the power W out carried by the flow and the heat production Wjj inside it: 

W m = W D + W out . (8) 

If the height of the granular column is not too small, we have, 

The contribution W[^ ct ' accounts for the kinetic energy of the fluctuating motion of the 
grains, which we assume is much smaller than the first contribution and may hence be 
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ignored; p ou t(< p) is the density of the orifice flow. For ideal liquid, Wd = 0, and the 
balance is given between Wi n and W out , resulting in the square root dependence Q ~ y/F of 
Torricelli. But for GOF, dissipation is so strong that the balance is mainly between Wi n and 
Wd- With the present experiment, (assuming p out = p/2,) we estimate W out /Wi n < 10 -2 , 
indicating that W ou t is negligible. In this situation, the rate Q (F) is determined by the 
dependence of the heat production power on the pushing force, Wd (F). Tailor expanding 
and comparing the result with (pQ), we have 

W in = W D = W D 1] F + wg ] F 2 + ... (10) 

and 

A = n ( p 2 pW D 1) /4g, (11) 

a = W D 2) /W D 1] . (12) 

Further study needs to consider the heat produced in the system. Employing Granular 
Solid Hydrodynamics (GSH) as given in [161 ]. we show a order-of- magnitude estimate for 
W D 1] below. 

First, note that in the present setup, dissipation occurs only in a small volume Vd close to 
the orifice at the bottom. Assuming that this region is a dense granular gas with a granular 
temperature T g , we have, according to GSH, the following relations for the pressure and 
heat production (up to quadratic order in T g ) 



ap pboTl 



2p c (l-p/ Pe ) 



l-a 



(13) 



W D = V D ioT g 2 (14) 

where p c ~ 1.0667p is the random closest packing density, the exponent a ~ 0.1, and &o ; 7o 
material parameters, of granular matter. Eliminating T g with Eqs.( ll3|[T4"]) and noting that 
P ~ AgF/TT(f) 2 we get 

W D = V D y/-^g^F (15) 
air(f) 2 p 



where 

pb (1 - gJo c ) a 



7 9 = 7 n 7 ° , ,a (16) 
pbn(l-oLPr) 



can be considered as a constant of about 0.5 Hz [l6j]. Inserting W D as given by (Tl5ll into 
dUD we get 

A = 2 lg0 V D (p c - p) /a. (17) 
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This A is comparable to the measurements if we take Vd as 10 times the orifice volume 
D 3 . To estimate a, higher order expansion terms are needed. The detailed analysis will be 
presented separately. 

Orifice flows are generally pushed by the force difference between the in- and outside 
of the orifice. The increase of the flow rate with the force obeys the Torricelli's law if an 
arch is absent, such as in the case of Newtonian liquids. Granular materials form FFA, the 
cause of the linear increase of flow rates measured in this work. The increase has not been 
observed before, perhaps because it is obscured by the Janssen effect. Our result shows that 
neither the minimum energy theory nor the hour glass theory is fully satisfactory, and the 
common view that the FFA is able to resist all loads on it is at most an approximation. 
As dissipative mechanisms are probably responsible for the behavior of dynamic arches, 
analyzing heat production should help to solve the long-term mystery of GOF. Moreover, 
we observed that pushing forces favorite a jamming transition, an effect also operative in 
the slight increase of the empty annulus size D a with the force. 
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